The role of the alloy structure in the magnetic behavior of granular systems 
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The effect of grain size, easy magnetization axis and anisotropy constant distributions in the 
irreversible magnetic behavior of granular alloys is considered. A simulated granular alloy is used to 
provide a realistic grain structure for the Monte Carlo simulation of the ZFC-FC curves. The effect 
of annealing and external field is also studied. The simulation curves are in good agreement with 
the FC and ZFC magnetization curves measured on melt spun Cu-Co ribbons. 

S: 

I. INTRODUCTION 

o 

Besides the importance of studying still open problems on basic magnetism, nanocrystalling systems attract more 
and more interest due to their applications on chemical catalysis and magnetic recording!! For the latter, the 
progressive increase of recording density has lead to the production of materials with smaller single domain particles. 
This requirement has the serious drawback that the effective magnetic moment of the particles suffer strong instabilities 
• of thermal origin, so called superparamagnetic limit. Also, when one deals with the nanometer scale, the magnetic 
systems are not easily reproduced and characterized, introducing additional difficulties for experimental studies. 
Experimental and theoretical results over the past years show that there are clearly many factors that can influence 
the magnetic and magnetotransport behavior of these systems, namely: The distribution of grain sizes, the average 
size and shape of the grains, the magnetic anisotropy of the individual grains and magnetic interactions among the 
nanometer-sized crystallites. 

The effect of some of these parameters on the magnetization and magnctorcsistancc has been partially investigated 
by several authors. El-Hilo et alu have used Monte Carlo simulations for determining the magnetoresistance depen- 
dence on the mean intergrariular distance, or rather, the particle concentration, using a simple expression previously 
obtained by Gittleman et alu. The influence of the log-normal distribution of magnetic moments and the ratio of the 
boundary toJhe volume scattering cross-sections on the magnetization and magnetoresistance has been examined by 
Ferrari et ala. Hickey et alu and Wiserlj have developed a phenomenological model which explains the almost linear 
variation of the magnetoresistance with the magnetization at low temperatures as a possible consequence of existing 
correlations between blocked and superparamagnetic particles. The blocking. effect on the. magnetization curves has 
been simulated by a Monte Carlo method in a paper of Dimitrov and WysirJj. Allia et alu have proposed analytical 
models that take explicitly into account the correlation arising from the dipolar interactions on nearly superparamag- 
netic systems. Also for interacting systems, Chantrcll et alB calculated the susceptibillity and Field-Cooling (FC) and 
Zero Field-Cooled (ZFC) curves for superparamagnetic particles, and Pike et alEl investigated the role of magnetic 
interactions on low temperature saturation remanence of fine magnetic particles. However, most of these models 
assume non-realistic size distributions of the grains, and constant anisotropy, as if all grains had the same shape. 
Both parameters originate effects that certainly are interesting to elucidate. 

In order to investigate the effect of different structural and magnetic factors on the magnetization properties of 
granular alloys we have simulated, by means of a standard Monte Carlo method, one of the most common character- 
ization techniques, the FC-ZFC curves. The samples were simulated by a Cell Dynamical System model which gives 
the realistic grain sizes and shapes distributions and, with the aid of experimental results, these quantities can be 
assigned to magnetization and anisotropy for each grain. After these quantities are known, the magnetic properties 
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of the system are studied as a function of temperature, applied field and thermal treatment of the sample. 

Comparisons between the realistic distribution and others are made to check the influence of grain sizes and shape 
distributions on ZFC-FC curves. An interesting correspondence is made between annealing and time iteration in 
CDS. The applied field effects on blocking temperatures and magnetization values for ZFC-FC curves are studied and 
agree with experiments. Our results are compared with experimental results also included is this work. 



II. SIMULATION OF THE GRANULAR ALLOY 



Granular materials correspond to metastable states resulting from the slowing down of the segregation process in 
binary off-critical mixturestil. In the case of metal alloys, a quench to the metastable region of the phase diagram 
generates a virtually permanent granular state at room temperature. Since the grain structure is not an equilibrium 
state, it cannot be derived from a minimization procedure. There are several statistically equivalent granular structures 
corresponding to different paths in phase space and one must go through the process of phase separation in order to 
reach one of these intermediate states and obtain realistic simulations of axanular alloys. 

The traditional approach to this problem is the Cahn-Hilliard equationu 

dj;(r,t) _ TT72 6Fty(r,t)] 
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where tp(r,t) is a conserved order parameter (in our case, the difference between Co and Cu concentrations), L 
is a phenomenological parameter related to the mobility, and F[ijj(r,t)] is the coarse-grained Landau free-energy 
functional: 



F[^(r,t)] = Jdr 



(2) 



with r and g positive phenomenological parameters. In principle, the solution of equation ([!]) yields the knowledge 
of the whole phase separation process. The problem is, of course, that ([[]) does not have an analytical solution and 
it is very hard to solve numerically. Almost all efforts in this sense were able to describe the very early .stages of the 
process only. Other common appro ach i es . ai j e My pte Ca rlo simulations with Kawasaki exchange dynamicsOEJ and Cell 
Dynamical System (CDS) modeling^E§lHllla'Ej'E2lEifL3. Since Monte Carlo simulations are computationally intensive 
and CDS has been successfully used to model this type of problems, we have chosen the latter. 

The basic point of CDS modeling is the discreteness of space and time. The point of view is equivalent to the 
one in the Cahn-Hilliard equation, in the sense that a coarse grained description based on densities is used, and 
the parameters are phenomenological. The stability of the dynamics and the computational efficiency allows the 
achievement of the asymptotic regime even for reasonably large systems. The disadvantage of this approach is that 
the arbitrary definition of parameters does not yield absolute information about the system, although this can be 
done in some cases by recovering the correct time and space scaleaB. 

The simulation, in this case, mimics an experiment in which a binary alloy is initially prepared at a temperature 
above the phase separation critical temperature, with given amounts of each component, and then quenched to a region 
of the phase diagram in which the homogeneous phase is metastable. The system is represented by a d-dimensional 
lattice and its configurations at time t and t + 1 are directly related by a map. Here we have considered only d = 2, 
since we believe that most of the properties we seek will be manifested in this dimension, so the computational effort 
involved in a 3d simulation would not improve significantly our results. 

We assume that the dynamics of each lattice cell is governed by a local relaxational racphaaism driven by a suitable 
map /. The exact form of the map is not important, as long as it has the correct flowEj£3o. For the segregation of 
binary mixtures we seek a map with one unstable fixed point at the origin and two hyperbolically stable fixed points 
at symmetrical positions. The stable fixed points correspond to the segregated phases rich in Co and Cu, and the 
unstable fixed point, to the homogeneous phase. 

The single cell dynamics is described by 

iP(t + l,n) = f[iP(t,n)], (3) 

where ijj(t, n) is the value of the order parameter in the cell n at time t. The addition of a diffusional coupling to its 
neighborhood leads to a non conservative dynamics of the form 

V(t + l,n) = fW,n)]+D[((tl>{t,n)))-fl>(t,n)] (4) 
= ^(t,n) + I{t,n) , (5) 
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where I(t,n) = f[i/j(t,n)] — i/j(n,t) + D [((ip(n,t))) — ip(n,t)] is the increment in order parameter after one iteration, 
and D is a positive parameter proportional to the phenomenological diffusion constant. ((*)) — * is the isotropic 
discrete Laplacian. We use the following definition of spatial average ((*)) on the two-dimensional square lattice: 

(m,n))) = ±£ ^ + ^E ^> ( 6 ) 

nn nnn 

where the sums are over nearest-neighbors and next-nearest-neighbors cells respectively. The conservation of order 
parameter may be locally imposed by the subtraction of the average gain in order parameter in the neighborhood of 
each cell. The conserved order parameter CDS model then reads 

i/>(t + 1, n) = V(t, n) + I(t, n) ~ ((I(t, n)» . (7) 



A. Simulations conditions 



We used a two-dimensional lattice with 450 2 cells and periodic boundary conditions. The order parameter represents 
the difference between the average concentration of Cu and Co in the alloy, ip = ipcu — V'Co- Our initial condition is 
the homogeneous system with ip(r) = tpo + S(f), where 5{r) is a small random fluctuation uniformly distributed in 
the interval [—0.005,0.005]. For a CugoCoio alloy, (V>Cu) = 0.9 and ("0Co) = 0.1 (average values of concentration), so 
we have ipo = 0.8. For such values of concentration the system is in a metastable state after the quench, so phase 
ordering proceeds via nucleation, and segregation of phases can occur only if nucleation centers are present. In our 
simulation nucleation centers consisting of four sites each with ip = —0.5, adding up to_about 10% of the lattice were 
randomly placed. For a longer discussion about off-critical quench in CDS models seeEEl. We used D = 0.5 and our 
choice of map is 

f(x) = .Atanhz (8) 

where A is a measure of the quench depth (A > 1 corresponds to a homogeneous mixture). In our simulations we 
use A — 1.3. Figure ^ shows the snapshots of the lattice after 100, 1000, 5000, 10000 iterations. The color black is 
associated to the presence of Co, so, according to our definition of "0, the grains have ip < 0. The time evolution of 
the pattern is driven by the gradient of the chemical potential, leading to domain coarsening and coalescence, so that 
patterns obtained after different number of iterations have distinct distributions of grain sizes and shapes, leading to 
distinct magnetic properties. As the segregation process proceeds, grains coalesce and grow, such that the number of 
grains decrease, as they become larger. In real samples such evolution of the grain structure is possible only with a 
thermal treatment. In this way, we are going to compare patterns simulated at later times with experimental samples 
that have been annealed at higher temperatures. 



B. Analysis of the grain structure 

We have used the standard algorithm developed by Hoshen and KopelmanEi to obtain the list of cells belonging to 
each grain. In order to check the form of the grain size distribution we have calculated the histograms corresponding 
to granular patterns after given number of iterations of Eq. [7] in terms of the magnetic moments m of the grains. As 
explained below, the magnetic moments are proportional to the number of sites in each cluster with -0 < 0. Figure 
H shows the histograms after 1000 (Sample 1), and 10000 (Sample 2) iterations, showing that as time goes on there 
is a shift to larger values of magnetization, due to the coalescence of grains. In order to compare the distribution 
obtained in simulated samples to the usually adopted log-normal distribution, we have fitted the above histograms 
with the curve 
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where P(m) is the probability density to find the grain magnetic moment between m and m + dm, m c is the most 
probable value of m, and w is the width of the distribution. 
The resulting fitting parameters are summarized in Table [I]. 
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FIG. 1: Segregation patterns for a Cu 90 Coi alloy after (a) 100, (b) 1000, (c) 5000, and (d) 10000 iterations of Eq. (@). Darker 
shades of gray correspond to Co rich regions. 
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FIG. 2: Histograms for the patterns after 1000, (a), and 10000, (b), iterations shown in Fig. |l| 
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sample 


m c {fJ>B) 


w 


A{(1 B ) 


1 


823 


0.80 


116 x 10 3 


2 


2091 


0.63 


189 x 10 3 



TABLE I: Parameters resulting from the fitting of a Log-normal distribution defined by Equation 
Fig- | 
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FIG. 3: Scatter plot of the anisotropy, measured as the ration Ri between the larger and the smaller dimension of the i-th 
grain, as a function of the number of sites, m in the grain. It is clear from this plot that the segregation process leads to 
larger and rounder grains, and that anisotropy is correlated to grain size. We expect the same behavior in real systems, that 
is, annealing will generate a system with larger but less anisotropic grains. 



Cobalt grains have fee symmetryEj, so it is reasonable to consider that the dominant contribution to the anisotropy 
energy comes from the grains shape. Observing the pattern formed after 1000 iterations (Fig. [l].(b)), for example, we 
notice essentially two kinds of grain shapes, one almost circular, and another elongated. Once the grains are labeled 
we can obtain information about their shape in order to calculate the direction of the easy magnetization axis. The 
method for calculating this direction is quite simple, it consists in calculating the position of the central site of the 
i-th grain, and then summing up the position vectors of each site related to the position of this central site. The 
resultant vector direction is the easy axis direction, e,. 

Finally we have translated our arbitrary units to physical ones, using a real CubgoCoio sample for which the average 
values of the grain magnetic moment and anisotropy constants are known. The determination of those experimental 



values is explained in Section II C. The first step was to associate a magnetization value to each grain. Assuming that 
each grain magnetization is proportional to its size, we can adjust our results by comparing them to experimental 
data. Of course one serious difficulty is that the whole simulation is 2d, so actually we have areas and not volumes 
for the grains. 

The average grain size is just the average number of lattice sites of the grains, n, and was associated to the average 
experimental magnetic moment m exp of a Cubgo Coio sample as 

nf m = m expo , (10) 

so that the magnetic moment of the i-th simulated grain is then m, = riif m , where m is the number of sites of that 
grain, and f m the conversion factor in units of m exp . 

We also need the conversion factor for the anisotropy constant K. We used the fact that, for cobalt, when R4 = 3.5 
the shape anisotropy constant is equal to K exp = 4.5 x 10 6 erg/cam 3 c3. With this it is possible to define the 
conversion factor fx as 

so that, for each grain, the anisotropy constant in the same units of K exp may be calculated as K{ — fuRi- 
The properties of the samples used in the Monte Carlo simulations are summarized in Table [n| 
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sample 


number of iterations 


N B 


m(fiB) 


A"(10 6 erg/cam 3 ) 


1 


1000 


539 


869 


1.58 


2 


10000 


371 


2140 


1.58 



TABLE II: Properties of the simulated CubgoCoio samples used in the Monte Carlo simulations. N g is the total number of 
grains in a sample, m(pfl) is the average magnetic moment of the grains and K is the average anisotropy constant of the 
sample. 



C. Experiment 

Continuous ribbons of CoioCubgo were obtained by melt spinning in AR atmosphere using a Cub-AR drum. The as 
cast material was subjected to furnace annealing at temperature range 400-600 °C for 60 min, generating samples with 
different panostructures. Details of the sample structure and how it is affected by the thermal treatment can be seen in 
referenceEZI. The magnetization measurements were performed on a commercial Quantum Design PPMS system with 
temperatures in the range 2-300 K, for different values of the applied field H. In each curve the sample was cooled 
to 2 K in zero field. After applying a magnetic field, the magnetization was measured with increasing temperature 
up to 300 K (ZFC magnetization). Subsequently, the magnetization was measured for decreasing temperature in the 
same field down to 2 K (FC magnetization). The magnetic field was applied by a superconducting magnet operating 
in persistent mode and the total time of each complete ZFC and FC curves was of 7 hours with a rate of temperature 
sweep of 1.5 K/min. 

The moment distribution of the samples was obtained by fitting a sum of Langevin functions weighted by a log- 
normal moment distribution in the m(H) curve measured at T = 300 K □. The saturation magnetization was obtained 
by extrapolation of the curve m(l/H) to zero temperature. 



III. MONTE CARLO SIMULATION 



A ferromagnetic particle becomes a monodomain when its linear size is below a critical value LL determined by the 
minimization of the total energy, including magnetostatic, exchange and anisotropy contributionsEj. Below this critical 
size, the energy associated to the creation of magnetic domain walls is larger than the decrease in the magnetostatic 
energy due to the smaller total magnetization. Such monodomain ferromagnetic particles can be viewed as large 
magnetic units, each having a magnetic moment of thousands of Bohr magnetons. Usually, in low concentration 
alloys, neighboring particles are separated by 10 — 30 nm, and direct exchange, as well as indirect (like RKKY) 
between particles is neglcctedcj. Thus, the magnetic properties of an assembly of nanoparticlcs are determined by 
the dipolar interaction energy between the particles along with thermal and magnetic anisotropy energies. In this 
paper we consider only the latter effect. Experiments conducted on magnetic nanoparticles show irreversible magnetic 
behavior below a irreversibility line Ti rr (H). In particular, the ZFC and FC magnetization curves do not coincide, 
and magnetic hysteresis appears. In such systems the origin of irreversibility is the interplay between thermal energy 
and some energy barrier, which hinders relaxation towards equilibrium. The magnetic irreversibility in nanoparticles 
is conventionally associated with the energy required—for a particle moment reorientation, overcoming a barrier due 
to shape, magnetoelasticity, or crystalline anisotropynS. 

Our system consists of N g magnetic monodomain particles, whose sizes, shapes and anisotropy constants are 
obtained by means of the previously described procedure. Each particle is described by its magnetic moment rftf, the 
direction of the easy magnetization axis, ej and its anisotropy constant Ki. All those quantities result from the CDS 
simulation of the granular structure. In the presence of an external magnetic field H , the Hamiltonian of our system 
can be written as 



-rft, ■ i? 



rfti ■ e. 



(12) 



where Ki = KiNi, Ni being the number of atoms of grain i. 

The ZFC-FC curves correspond to noncquilibrium states of the system, therefore they are dependent on the tem- 
perature variation rate. In terms of a Monte Carlo simulation this means that we cannot wait too long at each 
value of temperature, and the usual mechanism of time averaging instead of ensemble averaging is not valid. The 
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estimation of the number of Monte Carlo steps in each temperature was completely arbitrary and it has no obvious 
relation to the actual time used in experiments. However, as we keep this number constant in our simulations we can 
at least guarantee that the cooling and heating rates are equal in all simulations and our conclusions will be valid 
at these rates. The averaging was done over 200 independent but statistically equivalent sweeps. For the ZFC curve 
we started from a configuration where the magnetic moments of the grains were randomly chosen and H = 0. The 
simulation proceeded by turning on the external field (typically fO-fOO Oe) in the x— direction at low temperature 
(T = 2 K). The N g grains were then sequentially chosen, and had their magnetic moment rotated by an angle sorted 
from a uniform distribution. The change in energy (AE) was calculated and the rotation accepted with probability 
p = min[l,exp(— AE/ksT)]. The update of the N g particles was repeated 100 times for the initial value of temper- 
ature, and then the temperature was increased by 2 K, and so on. At T — 200 K we kept the last configuration of 
magnetic moments and used it as the initial configuration for the FC curve, following the same procedure used in the 
ZFC curve, only decreasing the temperature 2 K each step. 

IV. RESULTS 

The analysis of ZFC curves usually involves two temperatures Tm and T irr defined as the temperature at the 
maximum and the temperature above which the system shows thermodynamic equilibrium properties corresponding 
to a superparamagnetic behavior, respectivelyEl For T < Ti rr irreversibility impedes the coincidence of the ZFC and 
FC curves. For T > Ti rr , the relaxation time for the magnetization of the largest particle is much smaller than the 
typical measuring time, and the ZFC and FC curves coincide. Our goal is to understand how these temperatures are 
affected by the distribution of sizes and shapes of the particles, magnitude of the applied field and annealing of the 
grains. 

A. Influence of the distribution of grain sizes and shapes 

In order to compare the effect of the distribution of grain sizes and anisotropics on the ZFC-FC curves, we performed 
a series of simulations using different types of distributions namely, uniform, Gaussian and realistic (obtained from 
CDS simulation) for the sizes, anisotropy axis sorted from a uniform distribution of directions, and obtained from 
the CDS simulation, anisotropy constants equal to the average value, and calculated from the simulated alloy. One 
important point in the choice of the anisotropy constants Ki, is that they are intensive quantities. In order to keep 
the total anisotropy energy independent of the choice of the K distribution, we calculated the average value of the 
anisotropy constant using the grain sizes as weights. The results below show the combination of those possibilities for 
the same sample under an applied external field H = 0.1 KOe. All curves show the behavior of the reduced average 
magnetization per grain, m/m„, where m s is the saturation value of the sample. 

We begin our study by considering a sample with the same number of grains, average magnetic moment and average 
anisotropy constant of sample 1, as defined in Table |n[ but with arbitrarily chosen size distributions. 

First we assign a constant value of magnetization, rrii = m = 869/is, and uniform anisotropy constant, Kt = K = 
1.58 x 10 6 erg/cm 3 , to all grains and anisotropy axis sorted from a uniform distribution. These choices represent a 
system in which all grains have the same size and shape, with anisotropy axis determined by the random alignment of 
the crystalline axis and the magnetic field. Figure |^ shows the results for this system. We notice that T irr = Taj = 8 K 
and that the peak of the ZFC is sharp. 



Next we improve the description by considering a system with magnetizations sorted from a Gaussian distribution. 
In order to have a distribution quantitatively similar to the simulated one, we have chosen it to be centered at 
m = 869/is with variance corresponding to the width of the distribution of sample 1, that is, a = 550/1^. We still 
keep the constant value of anisotropy K, i.e. all particles with the same shape, and anisotropy axis sorted from the 
uniform distribution. With these choices we have a system with particles of different sizes, but uniform shape, and 
with easy magnetization axis determined by the grain crystalline axis, as before. As can be seen in Figure |5J, the 
existence of particles with different sizes is not a sufficient condition for having Ti rr ^ Tm as observed in experimental 
curves. A symmetrical distribution of sizes leads to curves similar to the ones obtained from a system of particles 
with uniform size in the sense that T irr = Tm in both cases. However, because of the inclusion of a size distribution, 
we have now grains larger than in the previous case, and this is reflected in the larger values of those temperatures. 
Bigger particles with higher magnetizations weight more in the ZFC curve than the smaller ones, and because of its 
bigger size they yield bigger values of T irr . Therefore, the position of T irr and Tm shift to higher temperatures and 
now Ti rr = Tm = 12 K. 
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FIG. 4: Simulated ZFC-FC curves for a sample of 539 grains with constant magnetic moment (= 869/j.b) and anisotropy 
constant (K = 1.58 erg/cm 3 ) and anisotropy axis sorted from a uniform distribution. Since all particles have the same size, 
T irr = Tm(= 8K in this case), as expected. 
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FIG. 5: Simulated ZFC-FC curves for a sample of 539 grains with Gaussian distribution of magnetic moments (m = 869 /is 
and a = 550 /is). All grains have the same anisotropy constant (K = 1.58 erg/cm 3 ) and the anisotropy axis are sorted from a 
uniform distribution. Although grains of different sizes are present, we still have Tj rr = Tm, indicating that the existence of a 
size distribution may not be the sole reason for having T rr > Tm- Here Ti rr = Tm = 12 K. 



The effect of the distribution profile can be clearly seen if we use the asymmetrical distribution generated by the 
CDS simulation, together with constant K and uniformly distributed anisotropy axis. Figure ^ shows the resulting 
curves for this case, where Ti rr > Tm (Ti rr = 18 K and Tm = 13 K) and both are larger than the values obtained 
in the previous case, with the symmetrical distribution. Moreover the maximum value of m/m s in the ZFC curve 
is lower as compared to the symmetrical distribution curve. In this case, the number of smaller particles is bigger, 
as compared to the monodispcrse case. Some small particles are becoming disordered at the temperature of the 
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FIG. 6: Simulated ZFC-FC curves for a sample of 539 grains with a realistic distribution of magnetic moments from the CDS 
simulation in which m = 869 (ib- As in Fig. H all grains have the same anisotropy constant (= 1.58 erg/cm 3 ) and the anisotropy 
axis are sorted from a uniform distribution. In this case Ti rr — 18 K and Tm = 13 K. 
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FIG. 7: Simulated ZFC-FC curves for sample 1 of Table |ll|. The introduction of a realistic distribution of anisotropy constants 
generates a ZFC curve more similar to an experimental one (see Fig. ^|) , as the ZFC and FC curves separate at a temperature 
well below T M . Here T irr = 36 K T M = 8 K. 



maximum magnetization of the sample, therefore, the magnetization peak in the ZFC curve is lower. Also for this 
reason, Tj rr is larger. 

Finally we use the complete information from the CDS simulation, that is, distribution of sizes, axis and anisotropy 
constants calculated from the simulated grains. In this case, Ti rr = 36 K and Tm — 8 K. 

Fr om now on, all the simulated curves use the structural properties of the simulated alloy, as explained in Section 

ETbI. 
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FIG. 8: (a) Simulated ZFC-FC curves for samples 1 and 2. (b) Experimental curves for a sample annealed for an hour at 450 °C 
and 550 °C Higher annealing temperatures facilitate the process of segregation, so we have a correspondence to a simulated 
sample with a larger number of iterations of the time evolution equation of the segregation process, Eq. (Q). 



B. Effect of annealing 

Even for non interacting grains, annealing is capable of drastically changing the ZFC-FC curves. As the sample is 
brought to a temperature closer to the phase separation critical temperature, diffusion is enhanced and fabrication 
defects are relaxed. The overall effect is an isotropic sample with larger and fewer grains. For this reason we have 
chosen an annealed sample as our starting reference system. As the thermal treatment is repeated, grains grow larger 
and coalescence further decreases the number of grains. It is possible to achieve the same result in the simulated alloy 
simply by iterating more times the time evolution rule defined by equation (0) . For this reason we assume that sample 
2, collected 9000 iterations after sample 1 was collected, represents an annealed sample, as compared to sample 1. 
Using the same conversion factors defined for sample 1 in s ample 2, we guarantee consistency in the values of rrii and 
Ki. 

As expected for a sample with larger grains, both Tm and Ti rr are larger, and the maximum of both curves are 
higher. While a direct comparison is difficult, we can see that there is a qualitative agreement with experimental 
results, as illustrated in Fig. ||.(b). Because of the existence of bigger grains in sample 2, the mean value of the 
magnetization increases; therefore Ti rr increases and also the low temperature limit of m/m s on the FC curve. 
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FIG. 9: Simulated ZFC-FC curves for sample 1 with different values of the applied field. Both Ti rr and Tm decrease with 
increasing field as can be seen in Table [II 
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FIG. 10: Experimental ZFC-FC curves for a sample annealed at 450 °C. Corresponding values of Tj rr and Tm are listed in 
Table |. 



Effect of external field 



Here we consider the effect of varying the magnitude of the applied magnetic field on sample 1. Figure || illustrates 
our results for fields in the range 50-1000 Oe. Simulated curves have the same qualitative behavior found in experi- 
ments, as depicted in Figure |l(]. That is, both ZFC and FC curves are dislocated to larger values of m/m s as the field 
increases, and Tm and Ti rr decrease. Table [II shows an estimation of the values of Ti rr and Tm for the simulated 
curves. 
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Simulation 


Experiment 


field (Oe) 


Tirr (-^) 


Tm(K) 


Ti rr (K) 


T M 


50 


36 


8 


39 


16 


100 


36 


8 


38 


15 


200 


32 


8 


35 


14 


500 


26 


4 


22 


11 


1000 


14 


2 


19 


9 



TABLE III: Values of T lrr and T M for the curves in Fig. § and [lo[ 



Both values of temperature reflect the competition between thermal excitation and the energy barrier between the 
two axial orientations. Ti rr is strongly influenced by the external field. As Ti rr is related with equilibrium, i.e, it is the 
temperature above which the system behaves like a superparamagnet, by increasing the external field, the relaxation 

time diminishes. For superparamagnetic systems, the relaxation time at temperature T is given by r = tq exp ^^y^ , 

where tq is a constant of the order of 10 -9 s and AE is the energy barrier that each particle has to overcome to 
minimize its energy. Because the energy barrier decreases under an external magnetic field, so do Tm and Ti rr , and 
the equilibrium state can be reached at lower temperatures. 



V. CONCLUSIONS 



There are several difficulties involved in the direct comparison of simulated and experimental curves, mostly because 
we are dealing with nonequilibrium states. The first issue is how to identify a simulated Co x Cui_ x sample with a 
real one. As explained above, the CDS simulation deals with densities, and the space-time scale is undefined. As a 
zero order approximation we can force the simulated alloy to have the same average value of magnetic moment, but 
the variance of the grain size distribution in the simulation cannot be adjusted to be exactly the same as in the real 
system. Besides all this, there is also the problem of relating Monte Carlo steps with cooling and heating rates of the 
ZFC and FC curves, as discussed above. In summary, even though we have used experimental parameters to keep 
our simulated curves as close as possible to the real ones, a direct comparison of them is too ambitious. However, 
qualitative comparisons to determine the effect of each of the parameters of our model are valid, and the importance 
of relating anisotropy, size and shape distributions is stated. 

The existence of a temperature Ti rr ^= Tm marking the irreversible portion of the ZFC curve is already well 
understood to be a consequence of a grain size distributioriEj. What we observed in this work is that the sole 
existence of size distribution is not sufficient to produce a simulated ZFC curve similar to the experimental one. 
For symmetrical size distribution of the magnetic particles, Tm coincide with Ti rr , within numerical accuracy as can 
be seen in Fig. |^. However, for non symmetrical particle size distributions, Ti rr > Tm as shown in Fig. || The 
introduction of a distribution of anisotropy constants is a delicate point. As can be seen in Fig. ^, shape and size are 
correlated, so one cannot independently sort the values of magnetic moments and anisotropy constants, even when 
using the correct distributions. This effect has been so far neglected in simulations of granular systems. The qualitative 
comparison of the ZFC-FC curves corresponding to samples that have evolved in time with experimental curves for 
thermally treated samples reinforces this idea. As the binary alloy relaxes towards its equilibrium configuration, fewer 
and larger grains are present, and only the larger ones show appreciable anisotropy as shown in the scatter plot of 
anisotropy as a function of particle size for samples obtained after times (Fig. ^) for samples 1 and 2 obtained after 
different number of iterations. 
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